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Some new results on nonperturbative renormalisation of quark bilinears in quenched QCD with Schrodinger 
Functional techniques are presented. Special emphasis is put on a study of the universality of the continuum limit 
for step scaling functions computed with different levels of O(a) improvement. 



1. Step scaling functions of quark bilinears 

Nonperturbative renormalisation is a key in- 
gredient in precision lattice QCD computations. 
Schrodinger Functional (SF) techniques enable us 
to determine nonperturbatively the renormalisa- 
tion group (RG) running of fundamental param- 
eters (coupling and quark masses) and compos- 
ite operators in QCD over a vast range of scales, 
thus providing a good control of systematics in 
the renormalisation procedure T. One impor- 
tant issue in this context is the extrapolation to 
the continuum limit (CL) required to remove the 
cutoff dependence in the RG running. Whenever 
this dependence is steep and cannot be tamed 
by appropriate O(a) improvement of the action 
and/or the relevant operators, the extrapolation 
may turn into a major source of uncertainty. 
Good control of the latter is also interesting in 
order to test universality of the CL for different 
lattice regularisations, along the lines of e.g. 
Here we report on an ongoing detailed investiga- 
tion of the CL extrapolation in the determination 
of the RG running of quark masses and quark bi- 
linears in quenched QCD with Wilson fermions. 

To describe the nonperturbative RG running 
of a bilinear Or = ipi^^ji where i,j are differ- 



ent quark flavours, we use a step scaling function 
(SSF) er r , defined as 



or (u) = exp ■ 
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where 7 is the anomalous dimension of the opera- 
tor, f3(g(/j,)) = dg(fi)/d(hifj,), and g is the renor- 
malised coupling. In order to compute a SSF in 
a SF scheme, the theory is put to live in a finite 
box of size L 3 xT with SF boundary conditions 1 . 
The renormalisation scheme is fixed by setting 
/i = 1/L and imposing in the chiral limit a renor- 
malisation condition of the form 3 : 
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where F is a suitable SF correlation function, and 
O is a factor (see below for its precise form) which 
cancels the multiplicative renormalisation of SF 
boundary fields. 

The numerical determination of the SSF at a 
given value of u = g 2 starts with the computation 
of the renormalisation constant on lattices of sizes 
L/a and 2L/a at fixed bare coupling; this allows 
to construct the finite-cutoff SSF Sr(w, a/L) = 
Zr(g 2 ,a/2L)/Zr(g 2 ,a/L). Then or is obtained 

1 We will always take T = L, vanishing boundary gauge 
fields and a phase = 0.5 in spatial boundary conditions. 
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by repeating the procedure at various values of a 
for fixed L, and extrapolating to the CL: 



ar(u) = lim Ylr(u,a/L) 

a/L-^O 



(3) 



g2( L ~i) =u 



We will concentrate on two bilinear opera- 
tors, namely P = 75 and Tok = i O aQk , witl 
a ok = |[70jTfc]- To compute renormalisation con- 
stants, Eq. is imposed with F and 2 set to 
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/l = -^6(^iN^[75]) 

for the pseudoscalar density and 
M^o) = -Q^{$i{x)vok4> 3 {x)O r Alk\) 
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for the tensor bilinear, where Oy[r] = 
fl6 S y z Ci(y)rCi( z ) is a SF boundary source. The 
resulting renormalisation constants and SSFs will 
bear the labels P and T, respectively. 

We stress that the computation of ct p and ct t 
suffices to determine the RG running of all quark 
bilinears in the CL. Recall also that by using axial 
Ward-Takahashi identities it is easy to show that 
cr P = cr^j , where a m is the SSF for quark masses. 

2. Cutoff effects in step scaling functions 

Cutoff effects in E P and E T can be analysed by 
considering the Symanzik expansion of the two- 
point functions entering their definitions. In the 
absence of 0(a) improvement both quantities are 
hence predicted to exhibit leading cutoff effects 
linear in a. To implement complete on-shell 0(a) 
improvement at the action level, it is enough to 
include the clover term proportional to c sw and 
the boundary counterterms proportional to ct and 
Ct 0]. 0(a) improvement for the operators in the 
chiral limit can be written schematically as 

P 1 = P ; T 7 fc - T ofc + a c T [doV k - d k V ] , (8) 

with Vfx — tpi^/^ipj. The coefficient c sw is known 
nonperturbatively in the whole range of values of 



go needed, while Ct and Ct are known in perturba- 
tion theory to g^ and g^ orders, respectively 
Gr, on the other hand, has been computed non- 
perturbatively only at a few values of go 6 . 

We have computed the SSFs E P and E T in 
quenched QCD at 14 different values of the renor- 
malised coupling, ranging from g (L^^) — 3.48 
to 3 2 (i~ 1 1 n ) = 0.8873, and four different values 
of the bare coupling constant in each case (cor- 
responding to L/a = 6, 8, 12, 16). The details are 
essentially the same as in The main novelty 
is that we have used two different regularisations, 
involving different levels of 0(a) improvement, 
the results of which will be referred to as [I] and 
[U], respectively: 

[I] Clover action with nonperturbative c sw , 
one-loop values for ct and dt- 

[U] Wilson action (c sw = 0), c t = 1 (tree level 
value), one- loop value for c t . 

S T has been always computed with an unim- 
proved T ok . 

As the effect of neglecting O(#o) terms in c t 
and ct can be expected to result in very small 
0(a) effects in the SSFs 0, the results for Y$ 
should approach the continuum limit as a 2 ; s[> U ', 
on the other hand, is expected to exhibit a linear 
behaviour in a. This is indeed observed in nu- 
merical results, as exemplified by Fig. ^ where 
the CL extrapolation for er P at two different val- 
ues of the renormalised coupling is shown. By fit- 
ting [I] and [U] results independently with various 
Ansatze, universality of the CL has been checked, 
up to deviations which never exceed 1.5 ct and can 
be regarded as statistical. Eventually, [I] and [U] 
results can be combined into a constrained fit to 
obtain a more precise determination of ct p , as ad- 
vocated in [7] for computations where 0(a) im- 
provement of operators is difficult to implement. 

S T is expected, on the other hand, to exhibit 
linear cutoff effects even when it is computed with 
an 0(a) improved action, due to the lack of oper- 
ator improvement. Furthermore, Ep' may show 
large 0(a) effects. This is indeed observed in pre- 
liminary results for this SSF (Fig.|2J. In this case, 
the two-action strategy reveals its usefulness, as 
assuming universality of the CL and performing a 



3 



0.975 
0.97 
0.965 

^ 0.86 
5 0.84 



f- 0.82 



i i i i i i i 



"13 



I i ll il l In: 



1 1 1 1 1 1 1 



_L 



0.05 0.1 0.15 
a/L 

Figure 1. CL extrapolation of E P for g 2 = 0.8873 
(top) and g 2 = 3.48 (bottom). Full (open) sym- 
bols correspond to [I] ([U]) data. Fits, exclud- 
ing L/a — 6 points, have been performed with 
Ansatze of the form a + p(a/L) 2 and a + p(a/L), 
respectively. 
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Figure 2. Same as Fig. Q]fbr £ T , showing a con- 
strained CL extrapolation. Linear behaviour in 
(a/L) has been assumed both for [I] and [U] data. 

constrained fit yields much more precise CL val- 
ues than would follow from [I] data alone. The 
corresponding result for a T is depicted in Fig. 

3. Outlook 

The results obtained so far demonstrate uni- 
versality of the CL result for the RG running 
of quark bilinears computed with actions with 
different levels of O(a) improvement. A better 
control of the CL extrapolation can be achieved 
in cases where nonperturbatively O(a) improved 
operators are not available by combining results 
from different actions. This technique has been 
used in our preliminary computation of the non- 




Figure 3. Preliminary results for function 
of u = g 2 from a constrained fit. The solid line is 
the one-loop result for cr T . 

perturbative RG running of tensor bilinears in 
quenched QCD. Pending tasks include compari- 
son with cutoff effects in perturbation theory, and 
the computation of the NLO anomalous dimen- 
sion of Tot;, required to match our results for cr T 
to conventional rcnormalisation schemes. 
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